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A plane-symmetric non-static cosmological model representing a bulk viscous fluid distribution has 
been obtained which is inhomogeneous and anisotropic and a particular case of which is gravitationally 
radiative. Without assuming any adhoc law, we obtain a cosmological constant as a decreasing function 
of time. The physical and geometric features of the models are also discussed. 

1. Introduction 

In recent years cosmological models exhibiting plane symmetry have attracted the attention 
of various authors. At the present state of evolution, the universe is spherically symmetric 
and the matter distribution in it is isotropic and homogeneous. But in its early stages of 
evolution, it could have not had a smoothed out picture. Close to the big bang singularity, 
neither the assumption of spherically symmetry nor of isotropy can be strictly valid. So, we 
consider plane symmetry, which is less restrictive than spherical symmetry and provides an 
avenue to study inhomogeneities. Inhomogeneous cosmological cosmological models play 
an important role in understanding some essential features of the universe such as the for¬ 
mation of galaxies during the early stages of evolution and process of homogenization. The 
early attempts at the construction of such models have done by Tolman 1 and Bondi 2 who 
considered spherically symmetric models. Inhomogeneous plane-symmetric models were 
considered by Taub 3 and later by Tomimura, 4 Szekeres, 5 Collins and Szafron, 6 Szafron 
and Collins. 7 Recently, Senovilla 8 obtained a new class of exact solutions of Einstein’s 
equation without big bang singularity, representing a cylindrically symmetric, inhomoge¬ 
neous cosmological model filled with perfect fluid which is smooth and regular everywhere 
satisfying energy and causality conditions. Later, Ruis and Senovilla 9 have separated out a 
fairly large class of singularity free models through a comprehensive study of general cylin- 
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drically symmetric metric with separable function of r and t as metric coefficients. Dadhich 
et al. 10 have established a link between the FRW model and the singularity free family by 
deducing the latter through a natural and simple inhomogenization and anisotropization of 
the former. Recently, Patel et al. 11 presented a general class of inhomogeneous cosmo¬ 
logical models filled with non-thermalized perfect fluid by assuming that the background 
spacetime admits two space-like commuting killing vectors and has separable metric coeffi¬ 
cients. Bali and Tyagi 12 obtained a plane-symmetric inhomogeneous cosmological models 
of perfect fluid distribution with electromagnetic field. Recently, Pradhan and Yadav 13 
have investigated a plane-symmetric inhomogeneous viscous fluid cosmological models 
with electromagnetic field. 

Models with a relic cosmological constant A have received considerable attention re¬ 
cently among researchers for various reasons (see Refs. 14-18 and references therein). Some 
of the recent discussions on the cosmological constant “problem” and on cosmology with 
a time-varying cosmological constant by Ratra and Peebles, 19 Dolgov 20-22 and Sahni and 
Starobinsky 23 point out that in the absence of any interaction with matter or radiation, the 
cosmological constant remains a “constant”, however, in the presence of interactions with 
matter or radiation, a solution of Einstein equations and the assumed equation of covari¬ 
ant conservation of stress-energy with a time-varying A can be found. For these solutions, 
conservation of energy requires decrease in the energy density of the vacuum component 
to be compensated by a corresponding increase in the energy density of matter or radiation. 
Earlier researchers on this topic, are contained in Zeldovich, 24 Weinberg 15 and Carroll, 
Press and Turner. 25 Recent observations by Perlmutter et al. 26 and Riess et al. 2 ‘ strongly 
favour a significant and positive A. Their finding arise from the study of more than 50 
type la supernovae with redshifts in the range 0.10 < 2 < 0.83 and suggest Friedmann 
models with negative pressure matter such as a cosmological constant, domain walls or 
cosmic strings (Vilenkin, 28 Garnavich et al. 29 Recently, Carmeli and Kuzmenko 30 have 
shown that the cosmological relativity theory (Behar and Carmeli 31 ) predicts the value 
A = 1.934 x 10 -35 s -2 for the cosmological constant. This value of A is in excellent agree¬ 
ment with the measurements recently obtained by the High-Z Supernova Team and Super¬ 
nova Cosmological Project (Garnavich et al; 29 Perlmutter et al.; 26 Riess et al.; 2 ' Schmidt 
et al. 32 ). The main conclusion of these works is that the expansion of the universe is accel¬ 
erating. 

Several ansdtz have been proposed in which the A term decays with time (see Refs. 
Gasperini, 33,34 Berman, 35 Freese et al., 18 Ozer and Taha, 18 Peebles and Ratra, 36 Chen 
and Hu, 3 ' Abdussattar and Viswakarma, 38 Gariel and Le Denmat, 39 Pradhan et al. 40 ). Of 
the special interest is the ansatz A oc S~ 2 (where S is the scale factor of the Robertson- 
Walker metric) by Chen and Wu, 37 which has been considered/modified by several authors 
( Abdel-Rahaman, 41 Carvalho et al., 18 Waga, 42 Silveira and Waga, 43 Vishwakarma. 44 ) 

In most treatments of cosmology, cosmic fluid is considered as perfect fluid. How¬ 
ever, bulk viscosity is expected to play an important role at certain stages of expanding 
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universe. 45-47 It has been shown that bulk viscosity leads to inflationary like solution, 48 
and acts like a negative energy field in an expanding universe. 49 Furthermore, There are 
several processes which are expected to give rise to viscous effects. These are the de¬ 
coupling of neutrinos during the radiation era and the decoupling of radiation and matter 
during the recombination era. Bulk viscosity is associated with the GUT phase transition 
and string creation. A number of authors have discussed cosmological solutions with bulk 
viscosity in various context. 50-53 

Roy and Narain 54 have obtained a plane-symmetric non-static cosmological models in 
presence of a perfect fluid. In this paper, we will investigate a plane-symmetric inhomoge¬ 
neous bulk viscous fluid cosmological models in the presence of a variable cosmological 
constant varying with time. In our previous paper 13 we have obtained a non-degenerate 
Petrov type-II solution. In the present paper we have derived a cosmological model, which 
has in general, as Petrov type-I solution and as a sub case it represents a gravitationally 
radiating Petrov type-II solution. The physical and geometric behaviour of the models will 
be discussed. 


2. The Metric and Field Equations 

We take the plane-symmetric spacetime considered by Roy and Narain 54 

ds 2 = dt 2 — dx 2 — B 2 dy 2 + C 2 dz 2 , (1) 

where the metric potentials B and C are functions of x and t. The energy momentum 
tensor in the presence of bulk stress has the form 

T ij = (p + p)viVj - pgij (2) 

and 

P = p-Zv* i (3) 

Here p, p, p, and £ are the energy density, isotropic pressure, effective pressure, bulk vis¬ 
cous coefficient respectively and v t is the flow vector satisfying the relation 

g lJ ViVj = 1 (4) 


The Einstein’s field equations are 

Rij - l;Rgij + Rgij = SnTij (5) 

where A is the cosmological constant. Eqs. (2) and (4) for the metric (1) lead to 

v 2 = v 3 = 0 (6) 


= 1 


From Eq. (4), we have 


( 7 ) 
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The field Eqs. (5) for the line element (1) lead to 


O u , -V 2 , -l -®44 , C44 B1C1 — B4C4 
-8t T[(p + p) Vl + p] = — + —-- 

0 _ C 44 — C\\ A 

-8ttp =-—-A 

_ B 44 - -Bn 

-8ttp =---A 

of / , 2 -1 -Bn , Cn B1C1 - B4C4 

-8tT[( p + p)u 4 - p] = -g- + --+ A 


-87r[(p + p)77iu 4 ] = 


B 14 , C14 


LVr 1 ± / -L J 

The suffixes 1 and 4 at the symbols B and C denote partial differential with respect to x 
and t respectively. 


3. Solutions of the Field Equations 

From Eqs. (9) and (10), we have 


Buv _ C uv 

~W = ~cT 


i. , i. 

U= iz{x + t), v = 


From Eqs. (8)-( 13), we have 


Buu T B vv C uu + C V 
4 B + 4 C 


B U C V + B V C U 
2BC 


Let us assume 


Buu B vv Cuu CVV 
4B + 4 C 


B = a(u)./3(v) 


C = f(u).g(v) 


From Eqs. (13), (16) and (17), we have 

(au/a) = (gv/g) = , 
(fu/f) (P V /P) l! 

where k\ is an arbitrary constant. Integrating Eq. (18), we get 

a = Mf kl 


g = N P kl 


( 20 ) 
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where M and N are constants of integration. From Eqs. (15), (16), (17), (19) and (20), we 
obtain 


Uu/S) 

(fuu/fu) 


(k\ + l) 2 ^p + ki(k'l - l)j 


ki{kl - l)$f- + {k\{kl - l) 2 - {k\ + 1) 2 }|| 

where hi is an arbitrary constant. Integrating Eq. (21), we get 

/ = (au + 
and 

/3 = (mv + n) v 

where a, 6 , m and n are arbitrary constants and 


= kl 


h = 


k'2 


k\ + 1 


fc 2 — 1 ’ kf + 1 ki — 1 

By suitable transformation metric (1) reduces to the form 

ds 2 = dT 2 - dX 2 - u 2kl>1 V 2u dY 2 - U 2 ^V 2klV dZ 2 

where 

U = ^(X + T) and V = ^(X -T) 

The effective pressure and the energy density for the model (25) are given by 

Aki^iv 


8 p=- 


8np = 


(T 2 - X 2 ) 

4 pv(k\ + h + 1) 


■A 


-A 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(T 2 - X 2 ) 

For the simplification of £, we assume that the fluid obeys an equation of state of the form 


P = IP , 


(29) 


where 7(0 < 7 < 1 ) is a constant. 

Thus, given £(i) we can solve the cosmological parameters. In most of the investigations in¬ 
volving bulk viscosity is assumed to be a simple power function of the energy density 55-57 

m = a,/ oof 

where £0 and k are constants. If k = 1, Eq. (30) may correspond to a radiative fluid. 
However, more realistic models 58 are based on n lying in the regime 0 < n < 2. 


3.1. Solutions for £ = £0 
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In this case we assume fc = 0 in Eq. (30). Eqs. (30) and (27) become £ = £o = constant 
and 


8irp = 


8 tt £ 0 


( T 2 - X 2 ) 


1 + 


(. k\ + l ) 2 Akipu 

(fc 2 + 1) J “ (T 2 - X 2 ) 


-A 


Using Eq. (29) and eliminating p(t) from Eqs. (28) and (31) we obtain 


2pv(ki + l ) 2 _ 47 rg 0 T (fci + l ) 2 

(T 2 — X 2 ) (T 2 -X 2 )i [ (k'l + 1) 


(31) 


(32) 


3.2. Solutions for (; = £o p 

In this case we assume fc = 1 in Eq. (30). Hence, Eqs. (30) and (27) become 


t = Co P 


(33) 


and 


87 r 



Co P 

(T 2 - X 2 )? 


(fci + i) 2 y 

(*? + !)/. 


Ak\pv 
(T 2 - X 2 ) 


+ A 


Using Eq. (29) and eliminating p(t) from Eqs. (28) and (34) we obtain 


1 + 7 


Co 

(T 2 — X 2 )i 



(fci + i) 2 y 
(fc? + i )). 


A = 


(34) 


4/rz/ 

(T 2 - X 2 ) 


fci + (fc 2 + fci + 1 ) 


Co 

(T 2 ~X 2 )i 


(1 + 


(fci + l)% Y 

(fcr + i) V. 


(35) 


We have observed from Eqs. (32) and (35) that the cosmological constant is a decreasing 
function of time and it approaches a small value as time progresses (i.e., the present epoch). 
Thus, with our approach, we obtain a physically relevant decay law for the cosmological 
constant unlike other investigators where adhoc laws were used to arrive at a mathematical 
expression for the decaying vacuum energy. 


4. Some Physical and Geometrical Features of the Models 


We shall now give the expressions for kinematical quantities and components of con¬ 
formal curvature tensor. With regard to the kinematical properties of the velocity vector Vi 
in the metric (25), a straight forward calculation leads to the following expressions for the 
expansion 6 and shear tensor cry of the fluid. 


022 


On = 


(T 2 -X 2 )i 
2 fci 


1 + 


3 (fc 2 + 1) 


(*h + l) 2 l 
(^ + 1 ) . 

T 2 


(T 2 - X 2 )i 


(T 2 - X 2 )i 


2 fc 2 — fci — 1 
3 (fc 2 + 1) 


- fc(fci ~ 1 ) 


(36) 

(37) 

(38) 
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1 nv(ki — l) 2 
_ 12 [ UV _ 

_ 1 \v(ki - l){u(ki + 1 ) - 1 } /i(fci - l){n(h + 1 ) - 1 }' 

G ( 1313 ) - g [ y2 Jp _ 

1 n u {ki — l) 2 
_ 12 [ UV _ 

_ 1 [ Mfci - l) 2 ' 

(2323) 6 [ UV 

_ 1 \n(ki - l){/x(fci + 1) - 1} v{ki - l){v{ki + 1 ) - 1 }' 
(2424) ~~ 8 U 2 V 2 

1 fu/(ki - l) 2 

+ 12 [ UV 

_ 1 - l){v{ki + 1 ) - 1 } n(ki - l){/i(fci + 1 ) - 1 }' 

°( 3434 ) - 8 [-y2-^-. 


(45) 

(46) 

(47) 

(48) 
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1 

12 


mK^i - l ) 2 


C, 


(1414) 


uv 


- l) 2 


6 


UV 


c. 


(1224) 


/i{k\ - l){/i( fc i + 1) - 1} v(ki - l){v(ki - 1) - 1} 


U 2 


V 2 


c, 


(1334) 


M(k i - l){M fc i + 1) - 1} v(ki ~ l)M fc i + 1) - 1} 


U 2 


V 2 


(49) 

(50) 

(51) 

(52) 


From Eqs. (45) - (52) it is observed that if k\ = 1 the spacetime is conformally flat. The 
spacetime is, in general, of Petrov-type I. However if k\ ^ 1 and either p = fc| 1 | or 
v = fc 1 _ 1 then the spacetime will be of Petrov-type II. For this case the reality conditions 
p > 0 and p > 0 imply that 


8 k 2 

(fci + l) 2 {k 2 + 1)(T 2 — X 2 ) 


< A< 


4fci 

(fc 2 + 1)(T 2 — X 2 ) 


(53) 


For p = fci 1 ) _ 1 , we have 7/ = ^ fcl +l^y+i) and in this case the geometry of the spacetime 
(25) takes the form 


2 fci - 2 m 

ds 2 = dT 2 - dX 2 - uTUxnv {kl+1)(k i + 1) dY 2 - {/T^T+TTy^i + D^ + D^ 2 


(54) 


The physical components Chijk take the form 


Cab = 


/ 

0 
0 
0 
0 

V o 


0 


0 


0 


£-»? 0 0 

0 £ + 77 0 

0 0 2 £ 

0 77 0 

77 0 0 


where 


and 




77 = 


(*1 - l ) 2 


(hi + 1 ) 2 (k 2 + 1) 
kf(k 1 -1) 


(fcl + l)(fc 2 + l) 2 J 


0 

0 

V 

0 

< + v 

0 


1 

uv 

1 

y 2 


0 
0 
0 
0 
0 


(55) 


(56) 


(57) 


From Eq. (55) we conclude that the model exhibits gravitational radiation. For large value 
of X it gives a type-two null spacetime representing an outgoing radiation field, although it 
will not satisfy the reality conditions at X = 00 . The radiative term in Cab is the quantity 


77 = 


fci(fci-i) 


4 |_(fci + l)(fc? + l) 2 J 


1 

w 


5. Discussion and Concluding Remarks 
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We have obtained a new class of inhomogeneous plane-symmetric cosmological models 
with a bulk viscous fluid as the source of matter. Generally, the models are expanding, 
shearing, non-rotating and Petrov-type I non-degenerate in which the flow vector is geode¬ 
tic. Under certain conditions the models will be of Petrov-type II. It is well known fact 
that the free gravitational field affects the flow of the fluid by inducing the shear in the flow 
line. 63 It is therefore interesting to study the cosmological models with given Petrov-types. 
The cosmological constant in all models given in Sec. 3.1 and 3.2 are decreasing function 
of time. The cosmological consequences of a decaying cosmological term has been dis¬ 
cussed lucidly 35 , 59 “ 62 with an adhoc assumption of a decay law. However, in our models, 
we recover such a law without assuming any a priori law for A. Thus our models are more 
general than those studied earlier. 
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